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[Abstract] A novel technique for developing propellant slosh damping requirements 
with respect to the stability characteristics of large flexible launch vehicles is presented. A 
numerical algorithm is devised which allows an automated software program to rapidly 
converge to pseudo-optimal solutions that minimize required propellant slosh damping for 
multiple tanks while maintaining constraints on the frequency response characteristics of a 
particular open-loop plant transfer function. An implementation of the algorithm using a 
high-order linear model of the Ares I plant dynamics considers all relevant dynamic 
interactions of flexible body modes, propellant slosh, and nozzle inertia effects. A 
high-resolution propellant damping requirements table is produced that can be used for 
baffle design. The method is demonstrated to provide exceptional speed and accuracy when 
compared with the alternative human-in-the-loop approach. 

I. Introduction 

The influence of sloshing propellant on the stability and control of large flexible booster rockets is a well known 
problem that has been extensively researched throughout various manned and unmanned launch vehicle 
programs. 1,2 ’ 3 Excitation of liquid motion at its fundamental frequency can rapidly lead to an instability that renders 
the control system ineffective or results in excessive forces and moments that may compromise structural integrity. 4 
At the minimum, an unstable closed-loop propellant mode will result in an attitude limit cycle that can overstress 
TVC subsystems or unnecessarily expend maneuvering propellant. 5,6 It is of paramount importance that the liquid 
motion be well understood so that it can be appropriately stabilized or suppressed via active or passive techniques. 7 

A typical approach of the flight controls discipline is the use of a planar linear perturbation dynamics model of 
the vehicle that includes the coupled dynamics of the elastic modes and the sloshing propellant. 8,9 The external 
forces due to nonlinear propellant sloshing of each j of n tanks is approximated by a linear spring-mass-damper 
confined to motion normal to the vehicle thrust axis. The spring-mass-damper model is appropriate for rockets with 
cylindrical tanks where the dominant acceleration is normal to the liquid's lateral degrees of freedom. Other 
sloshing models, such as the pendulum approximation, may be necessary for certain vehicle or container 
configurations. 

The 2 nd -order linear slosh dynamics are described by the equations 

k 

+2 ^sj w sj b SJ + w)j 6 S] =- S-I sj (p - X <p siJ i), , j= 1 ...« (1) 

i= 1 

where 5 sJ is the perturbation state of an equivalent mass rn sj that approximates the portion of liquid that is 
involved in sloshing, cp is the perturbation angle relative to the trajectory, and 8- is the body lateral acceleration. 
The eigenvectors W Slj and the generalized accelerations h } represent forcing of the slosh dynamics via elastic 
motion. The non-sloshing fluid is assumed to be stationary in the body frame, and the equivalent natural frequencies 
co SJ and viscous damping % SJ are nonlinear functions of the axial acceleration and the liquid and container 
properties. 3 It is clear that the liquid motion is excited by rigid body translation and rotation as well as forces 
exerted by the elastic deflection of the vehicle at the tank nodes. The liquid response also drives both rigid body and 
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flexible modes. It is assumed that the perturbations are small and that the forcing function frequency content is 
sparse near the true natural frequency of the liquid motion. 2 3 

While it is possible to compute analytical stability criteria for reduced models of the launch vehicle dynamics, 6 
as multiple slosh and additional bending modes are included, the order of the model increases sufficiently to 
preclude an analytical solution. A numerical analysis of the frequency response characteristics of the high-order 
open-loop model can be utilized to assess the closed-loop stability with many of the relevant elements (actuator 
models, bending filters, and so on) included. Stability is characterized in terms of gain and phase margins in 
addition to other criteria that constrain open-loop response of the bending or sloshing modes. 8 

Sloshing propellant response is often neglected explicitly in the design of a controller and is passively 
suppressed using baffles or other mechanical devices. A necessary tradeoff is the sacrifice of injected mass to 
account for the increased mass of the baffle mechanisms versus additional control system complexity and risk 
associated with insufficient slosh damping. In the early design phases of the vehicle, the exact propellant damping 
characteristics and their interaction with the control system are largely unknown to the control designer. A baseline, 
conservative damping requirement may be specified, but must be iteratively modified as the vehicle design matures. 
The damping requirements at a particular liquid level may increase or decrease due to configuration changes or 
performance constraints, or the structure may not be able to accommodate a mechanical baffle at the specified liquid 
level. The ability to rapidly generate a damping requirements profile based on a fixed control system configuration 
can greatly accelerate the design lifecycle. 


II. Motivation 

Propellant slosh damping in the Ares I Crew Launch Vehicle (CLV) is of particular interest due to the potential 
for significant control-slosh interaction. The sloshing fluid is located far forward of the composite first stage center 
of mass, the vehicle is highly flexible, the axial acceleration profile varies greatly, and the vehicle and control 
system characteristics before and after staging are significantly different. The elastic deformation of the integrated 
stack during first stage ascent may couple with the liquid motion due to coalescence of the oxidizer slosh and first 
lateral bending mode frequencies at peak axial acceleration. Although the tank is nearly full and the liquid level is 
constant during first stage ascent, appropriate baffle placement is required to ensure gain stabilization of the oxidizer 
mode. An example of the time- varying nature of the slosh parameters is shown in Figure 1 . 




Figure 1: Typical slosh parameter variation 

During upper stage flight, the shift in frequency and increase in sloshing propellant mass as the tanks empty 
preclude gain stabilization of the slosh dynamics. Both slosh modes appear at frequencies near the rigid-body 
control crossover frequency, and while phase stabilization of the slosh response is straightforward, it naturally leads 
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to reduced phase margin if the response is underdamped. Most importantly, many of the system parameters have 
significant uncertainty. 

Accurate determination of acceptable propellant slosh response as it relates to vehicle stability margins and 
performance objectives is a difficult task. An optimal damping profile varies with liquid level as the liquid and 
vehicle mass properties change, and is further complicated by the presence of two interacting tanks whose 
cross-coupling effects may be significant. 

The compensated open-loop dynamics of the Ares I vehicle and similar configurations exhibit conditional 
stability at least in one dimension, and typically in both magnitude and phase. 8 As such, a stability analysis can be 
readily visualized using a Nichols chart, where the log magnitude of the proper plant open-loop transfer function is 
plotted with respect to phase on a set of orthogonal axes. A natural choice of loop opening occurs at the actuator 
input command, where the system dynamics are SISO and the classical gain and phase margins relate directly to 
multiplicative input uncertainty and delay uncertainty in the autopilot loop. An example of a typical Nichols 
response is shown in Figure 2. 


FRACTAL Ares I Stability Result 





pitch 


yaw 



phase 


Figure 2: Nichols response example 

The path through the resultant magnitude-phase space is traced in increasing frequency. The Nichols criteria 
are a direct mapping of the Nyquist criterion: knowledge of the location of the open loop poles and the number of 
encirclements of the critical point -1+ jO can be used to assess relative stability of the closed-loop system 
dynamics. Multiple gain and phase margins are computed, each corresponding to the approximate frequency of the 
associated physical dynamics. Other robust frequency- domain stability criteria may be applied to the open-loop 
frequency response, including a limitation on the peak open-loop bending mode response, shown as “flex peak” 
above. 

One approach to analyzing the stability of the closed-loop system utilizes the aforementioned baseline tank 
configuration. The robustness of the control system is determined by modifying system parameters either manually 
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or via a Monte Carlo method. If it can be determined that stability objectives are not achieved either in the time or 
frequency domain, the damping is increased in the time vicinity of the problem area and the simulation is re-run. 

This human-in-the-loop approach to the determination of propellant slosh damping requirements is slow and 
imprecise. A slightly more systematic approach allows an analyst to compare the frequency- domain trajectories of 
multiple systems with slightly varied damping characteristics. The analyst must select the best system response at 
the displayed time-slice. An example of this type of operator display appears in Figure 3. Clearly, the resolution is 
limited and the method is labor- intensive. 




phase 


phase 


phase 
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Figure 3: Manual damping determination method 

Any automated search of the solution space must attempt to determine, within reasonable bounds and with 
sufficient fidelity, the correct values of damping ratios for each tank at each time in the trajectory that result in 
acceptable time and frequency domain control system performance. It has been determined via simulation that the 
solution region is only vaguely known, that multiple minima exist, and that the solution is highly sensitive within a 
very small region. The computational requirements are daunting if the solution is to be obtained via a direct search 
with good precision. 

An efficient algorithm can capitalize on the ability of modem computer hardware to calculate the response of 
high-order linear dynamics models with exceptional speed as compared with computer hardware available for 
similar analyses in past launch vehicle programs. All dynamic characteristics can be considered, including multiple 
tanks and dozens of flexible modes, with minimal impact on computation rate. It is therefore possible to utilize 
high-order linear models within a specialized parameterization mechanism that can be used to automatically 
optimize various system parameters with greater speed and accuracy than possible with a human operator. The use 
of this capability in the design of a novel algorithm for optimizing propellant slosh damping requirements is the 
primary motivation of our approach. 


III. Problem Solution 

The proposed method of determining a near-optimal parameterization of the linear, time-invariant SISO system, 

x=A(0)x+bu , y=Cx (2) 

where 0 is a vector parameter, relies on a geometric interpretation of the complex frequency response that can be 
calculated as 

G(jco,e)=C[jwI-A(0)Y i b. (3) 

Let there be a function v that maps the complex valued function G(j to, 0) onto the Nichols plane such that 
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v(G{jw,0)) = 


j£.G(jw,0) 
201og|G( j to, 0 )| 


= rj (ju>,0) 


(4) 


which is the function of real frequency co appearing in a Nichols plot. Suppose that the gradient of r] with respect 
to the parameters 0 is known approximately or can be estimated, that is, the directions of the co -trajectory 
perturbations with respect to the k parameters of interest 0 Z can be determined: 




.... 


(5) 


Intuitively, the functions </>, point in the direction of the expected trajectory perturbation to be taken by the 
function r](jco,0) with respect to a small perturbation in . The values 4> , may be resolved, for example, into 
their magnitude and phase components along the frequency path. Intuitively, one may then prescribe convex 
regions S f m the Nichols plane whose boundaries have components that are nonnal to V 0 /] . Thereafter, define a 
binary function 


T(n,s)=\ l ^ s >’ v/ 

[ 0, otherwise 


( 6 ) 


which is presumably a binary test of the presence or absence of the intersection of the open-loop response r) subject 
to the parameterization 0 with the set of exclusion regions S . 

The binary function t , subject only to perturbations of the parameters 0 , is a straightforward test that may be 
incorporated into the following binary search algorithm, given initial guesses O imax [O] f 0 imin [ 0], 7=1 ... A: on the 
upper and lower bounds of each of the k parameters. The algorithm is formulated as follows: 


do [ n ] 

for each parameter 0, 
loop [ m ]: 
compute ri(0) 
if t(jj, S) 

i]=0fM 

0 ,.[ ot + 1 ] = 0 i .[ ot ]-( 0 ,.[/ w ]- 0 i . jBI < b [ ot ])/2 

otherwise 

d i [m+l}=0 i [m]+(0 imax [m]-0 i [m])l2 

end 

if 0,, m Jw + l]-0,, m/ „[m+l]<e,. 

0 i =O i [m\ + e i 

exit loop [ m ] 
end 

end loop [ m ] 
end for each parameter 0, 

f) = k O 7 = 1 k 

w i , min rv l w i,min> 1 

max ~~ ^2^i,max > ? _ 1 •••k 

while O imax -0 imi >e if i=\...k 

Pseudocode: Parametric binary search algorithm 


The algorithm progressively bounds the value of each parameter 0 ? so as to constrain the solution set of 
parameters 0 (where t= 1 ) within the arbitrarily small region e near the boundaries of S t . The algorithm is 
designed to converge efficiently to the semi-optimal values 0 7 when only one unique region exists with respect to 
0 j at the switching boundary of t , and it is known that the set of trajectory perturbation vectors (pji jcv) are 
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approximately aligned at similar co . Proper selection of the exclusion regions S f , the initial bounds 
1, 0 imax [0] , and the convergence criteria e, govern the algorithm's efficiency. 

A practical example applies the algorithm to a slosh response optimization problem. For a phase-stabilized 
slosh mode, it is postulated that the absolute magnitude of the open-loop gain shall not exceed some gain value 
g max . Therefore, the boundary of the region S { becomes the g max line in the phase-magnitude plane, and the 
remaining boundaries of the region are placed such that phase uncertainty is tolerated. An example of the phase- 
stabilized and gain-suppressed slosh exclusion region is shown in Figure 4. 



IV. Results 

The proposed algorithm is employed to determine a damping profile where the varied parameters 0 X 0 2 are the 
LOX and LH2 damping ratios, respectively. The system is most sensitive to LOX damping due to the LOX mode 
frequency and the large slosh mass, but both tank profiles can be optimized such that stability is maintained with 
minimum baffle mass. The use of the present technique allows rapid convergence to a near-exact solution. 
Typically, e ; = 2E-5 for an attained resolution in the damping ratio of approximately IE-4 . This is more than 
sufficient to specify for the purposes of baffle design. 

Algorithm convergence typically occurs in only 20-30 iterations and occupies about 30-60 seconds of computer 
time per trajectory time point on a modem (ca. 2008) desktop computer. The rapid computation cycle allows the 
determination of the optimum damping ratios with high resolution m the trajectory, typically every second. The use 
of a full-order linear model is employed so that all relevant flexible mode coupling effects are considered. 

While the region S x in Figure 4 is arbitrarily defined, it is intuitive that regions may also be defined with 
respect to stability margin objectives, so that the optimized response does not pass in the vicinity of a minimum gain 
or phase margin. The region S 2 shown in Figure 4 is such an example. 

Several points along the trajectory are considered, and the solution is determined at each point. The solution 
set, applied to a series of linear systems representing the open-loop controller-plant dynamics, yields a series of 
frequency response curves that ideally do not intersect the exclusion regions S , . An example of a solution set for 
upper-stage flight is shown in Figure 6. 
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Figure 5: Upper-stage trajectory solution set 

After suitable postprocessing, a minimum damping requirements table can be generated that catalogues the 
minimum acceptable damping values 0 L 0 2 that achieve the prescribed response characteristics with respect to the 
time evolution of the trajectory. This data is delivered to structures disciplines for use in baffle design. The time 
axis of the trajectory may be directly transcribed to nominal propellant and oxidizer liquid levels and used to 
generate a requirements table with respect to the fluid level in each tank. A tank liquid level versus damping 
requirements table is graphically depicted in Figure 6. 

Other extensions of the proposed algorithm have been implemented, including optimization of parameters in the 
presence of other modeling uncertainties. The approach relies on the speed of computing the linear system 
dynamics to systematically determine the required parameters that meet stability margin requirements while 
dispersing the remainder of the system parameters using Monte Carlo methods. Since the geometric criteria is 
binary, the algorithm need not compute the full set of random draws after the failure of a single test case, greatly 
reducing the overhead associated with the Monte Carlo method. Finally, the number of Monte-Carlo runs per test 
case can be increased using the relation such as 

w=«o+lylkr'] (7) 

where y is a gain parameter and n is bounded by n max , so that the full number of dispersions is not exercised 
until the solution is very near the converged value. 
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LOX Minimum Damping Requirement Profile 



damping requirement 


Figure 6: Minimum damping requirements example 


V. Conclusions 

The proposed method has been demonstrated to be effective in rapidly producing a set of minimum slosh 
damping requirements that satisfy flight control system performance objectives. The binary search algorithm 
capitalizes on known magnitude-phase response trends in order to construct pass- fail criteria and allow a solution to 
be determined much faster than with a direct search or a manual method. This technique has been used in other 
other parametric optimization and senstivity studies and will likely continue to reduce labor expenditure and 
enhance fidelity in the design and analysis of the Ares I and similar series of vehicles. 
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